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We present a selective overview of the current state of our knowledge (more pre- 
cisely of our ignorance) regarding the derivation of Fourier's Law, J(r) = — KVT(r); 
J the heat flux, T the temperature and k, the heat conductivity. This law is 
empirically well tested for both fluids and crystals, when the temperature varies 
slowly on the microscopic scale, with k an intrinsic property which depends only 
on the system's equilibrium parameters, such as the local temperature and density. 
There is however at present no rigorous mathematical derivation of Fourier's law 
and ipso facto of Kubo's formula for k, involving integrals over equilibrium time 
correlations, for any system (or model) with a deterministic, e.g. Hamiltonian, 
microscopic evolution. 

1 Introduction 

There are at least two distinct situations in which Fourier's Law is observed to hold 
with high precision: 

1. An isolated macroscopic system which is prepared at some initial time, say 
t = 0, with a nonuniform temperature To(r), e.g. a fluid or solid in a domain 
A surrounded by effectively adiabatic walls. At t > 0, the temperature will 
change, due to the heat, i.e. energy, current, with the energy density satisfying 
the conservation equation: 



where CviT) is the specific heat per unit volume and we have assumed that 
there is no mass flow or other mode of energy transport beside heat conduction 
(we also ignore for simplicity any variations in density or pressure). Eq.(|^) is 
to be solved subject to the initial condition T(r, 0) = To(r) and no heat flux 
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across the boundary of A. The stationary state, achieved as t — > c», is then one 
of uniform temperature T determined by the constancy of the total energy. For 
our purposes we can also think of A as a torus, i.e. having periodic boundary 
conditions. 

2. We consider the system in contact with heat reservoirs which specify a time 
invariant temperature Ta at points of the boundary r g (9A)q in contact with 
the a-th heat reservoir, a > 1. When the system has come to a stationary 
state (again assuming no matter flow) its temperature will be given by the 
solution of Eq.(|l]) with the left side set equal to zero, 

V • J(r) = V • (kVT (r)) = 0, (2) 

subject to the boundary condition T(r) — Ta for r e ((?A)q and no flux 
across the rest of the boundary which is insulating or periodic in the direction 
perpendicular to the heat flow. A simple example of this situation is the usual 
set up for a Benard experiment in which the top and bottom of a fluid in a 
cylindrical slab of height h and cross sectional area A are kept at different 
temperatures and T;, respectively. (To avoid convection one has to make 
Th > Tb or keep jT^ — Ti,| small). Assuming uniformity in the direction parallel 
to the vertical x-axis one has in the stationary state a temperature profile T{x) 
with f (0) = Tfc, f{h) = Th and k(T)^ =Const. for x e (0, h). 

From a physical point of view, which is how we presented them, the two cases 
are conceptually very similar (some physicists would even say identical). We have 
implicitly assumed that the system is described fully by specifying its temperature 
r(r, t) everywhere in A. What this means on the microscopic level is that we 
imagine the system to be in local thermal equilibrium (LTE). To make this a bit 
more precise we might think of the system as being divided up (mentally) into many 
little cubes, each big enough to contain very many atoms yet small enough on the 
macroscopic scale to be accurately described, at a specified time t, as a system in 
equilibrium at temperature T{ri,t), where is the center of the i-th cube. For 
slow variation in space and time we can then use a continuous description T(r,t). 

This notion is made precise in the so called hydrodpiajnic scaling limit (HSL) 
where the ratio of micro to macro scale goes to zerccZlLBEZl. The macroscopic 
coordinates r and t are related to the microscopic ones q and r, by r = eq and 
t = e"T, i.e. if A is a cube of macroscopic sides I, then its sides, now measured in 
microscopic length units, are of length L — e^^l. We then suppose that at t = Q 
our system of = pL'^ particles with Hamiltonian 
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is described by an equilibrium Gibbs measure with a temperature T(r) ~ T{eq): 
roughly speaking the phase space ensemble density has the form. 
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where Q = (qi, . . . ,qAr) € A'*^, P ~ (pi, . . . ,pAr) S M , 0(q) is sonie|Short range 
inter particle potential, u{(\i) an external potential and /3(7^(r) = ro(r)t^. 

In the limit e 0, p fixed, the system at i = will be macroscopically in LTE 
with a local temperature To(r) (as already noted we suppress here the variation 
in the particle density n(y)). We are interested in the behavior of a macroscopic 
system, for which e << 1, at macroscopic times t > 0, corresponding to microscopic 
times T = e~"i, a = 2 for heat conduction or other diffusive behavior. The implicit 
assumption then made in the macroscopic description given earlier is that since the 
variations in Tq{y) are of order e on a microscopic scale, then for e << 1, the system 
will, also at time t, be in a state very close to LTE with a temperature r(r, t) that 
evolves in time according to Fourier's law, Eq.(|l]). 

^From a mathematical point of view the difhcult problem is proving that the 
system stays in LTE for i > when the dynamics are given by a Hamiltonian time 
evolution. This requires proving that the macroscopic system has some very strong 
ergodic properties, e.g. that the only time invariant measures locally absolutely 
continuous w.r.t Lebesgue mfiasmie are, for infinitely extended spatially uniform 
systems, of the Gibbs type EEIIIj'Ej. This has only been proven so far for systems 
evolving via stochastic dynamics, e.g. interacting Brownian particles or lattice gases. 
In these systems the relevant conserved quantity is usually the particle density 
rather than the energy density. We shall not discuss such stochastic evolutions 
here but refer the reader to ||7^,0 for a mathematical exposition. 

The only Hamiltonian system for which a macroscopic transport law has been 
derived is a gas of nonintcracting particles moving among a fixed array of peri- 
odic convex scatterers (periodic Lorentz gas or Sinai billiard). For this system one 
can prove a diffusion equation like Eq.(|l]) for the density of the particles, both 
for the initial value and the suitably defined stationary state problem , .with the 
(self)diffusion constant given by the Einstein-Green-Kubo formula t2lEj. Unfor- 
tunately, the absence of interactions between particles makes this system a poor 
model for heat conduction in realistic systems. In particular there is no mechanism 
for achieving LTE. The speed of each particle |v| does not change in the course of 
time and the diffusion constant for each particle is proportional to its speed. The 
diffusion equation for the density mentioned above are therefore in fact separate 
uncoupled equations for particles with specified speeds. It corresponds to the usual 
diffusion equation only when all the particles have the same speed. 

To remedy this problem it would be necessary to add interactions between the 
moving particles, e.g. instead of points make them little balls, and then derive 
coupled equations for the diffusion of both particle and energy densities. This is 
what we would consider a satisfactory answer to the challenge in the title of this 
article and we offer a bottle of very good wine to anyone who provides it. We 
believe that this system, with only two conservation laws and an external source 
(the fixed convex scatterers) for chaotic dynamics may be the simplest Hamiltonian 
system for which such results could be proven rigorously. 

Just how far we are from such results will become clear as we describe our cur- 
rent mathematical understanding of the stationary nonequilibrium state (SNS) of 
macroscopic systems whose ends are, as in the example of the Benard problem, 
kept at fixed temperatures Ti and T2. The heat conductivity in this situation can 
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be defined precisely without invoking LTE. To do this we let J be the expecta- 
tion value in the SNS, i.e. we assume that the SNS is described by a phase-space 
measure (whose existence we discuss later), of the energy or heat current flowing 
from reservoir 1 to reservoir 2. We then define the conductivity as J /{A5T/L) 
where 5T / L = (Ti — Ti)lL is the effective temperature gradient for a cylinder of 
microscopic length L and uniform prossection A and k.{T) as the limit of when 
ST — > (Ti = T2 = T) and L — > oolj. The existence of such a limit with k positive 
and finite is what one would like to prove. 



2 Heat Conduction in Gases 



Before going on to a mathematical discussion of heat conducting SNS, we turn 
briefly to the "kinetic theory" analysis of heat conduction in gases. This is histori- 
cally the first example of a microscopic description of this macroscopic phenomenon. 
It goes back to the works of Clausius, Maxwell and Boltzmann B who obtained a 
theoretical expression for the heat conductivity of g Vt, independent of 

the gas density. This agrees with experiment (when the density is not too high) 
and was a major early achievement of the atomic theory of matteiEl. 

Clausius and Maxwell used the concept of a "mean free path" A: the average 
distance a particle (atom or molecule) travels between collisions in a gas with par- 
ticle density p. Straightforward analysis gives A ~ l/pira^, cr ^ an "effective" hard 
core diameter of a particle. They considered a gas with temperature gradient in 
the x-direction and assumed that the gas is (approximately) in local equilibrium 
with density p and temperature T(x). Between collisions a particle moves a dis- 
tance A carrying a kinetic energy proportional to T[x) from x io x + A/\/3, while 
in the opposite direction the amount carried is proportional to T{x -\- A-\/3). Taking 
into account the fact that the speed is proportional to a/T the amount of energy 
transported per unit area and time across a plane perpendicular to the x-axis J is 
approximately, 

J ^ py/r T(x)-T(x + \Vl) --(j-^VT—, (5) 
L J dx 

and so K ^ \/T independent of p, in agreement with experiment. 

It was clear to the founding fathers that starting with a local equilibrium situ- 
ation (corresponding to a Maxwellian distribution of velocities) there will develop, 
as time goes on, a deviation from LTE. They reasoned however that this deviation 
from local equilibrium will be small when (X/T)dT / dx « 1, the regime in which 
Fourier's law is expected to hold, and the above calculation should yield, up to some 
factor of order unity, the right heat conductivity. In fact if one computes the heat 
flux at a point x by averaging the microscopic energy current at cc j = pv(|mw^) 
over the one particle distribution function f{r,-v,t) then it is only the deviation 
from local equilibrium which makes a contribution. The result however is essen- 
tially the same as Eq.(^. This was shown by Boltzmann who derived an accurate 
formula for k in gases by using the Boltzmann equation to compute k. If one takes 
K from experiment the above analysis yields a value for a, the effective size of an 
atom or molecule, which turns out to be close to other determinations of the char- 
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acteristic size of an atomEI. This gave evidence for the reaUty of atoms and the 
molecular theory of heat. 

Using ideas of hydrodynamical space and time scaling described earlier it is pos- 
sible to derive a controlled expansion for the solution of the stationary Boltzmann 
equation describing thp steady state of a gas coupled to temperature reservoirs 
at the top and bottorrOOa. The coupling is implemented by the imposition of 
"Maxwell boundary conditions": when a particle hits the left (right) wall it get 
reflected with a distribution of velocities 



2T:{kTa 



\Vx\exp 



dw 



a = 1,2 



(6) 



corresf 
showsE 



g to a temperature Ti {T2) at the left (right) wall. One then 
^ that for e << 1, e being now the ratio A/L, the Boltzmann equa- 
tion for / in the slab has a time independent solution which is close to a local 
Maxwellian, corresponding to LTE, (apart from boundary layer terms) with a local 
temperature and density given by the solution of the Navier-Stokes equations which 
incorporates Fourier's law as expressed in Eq.(^. The main mathematical problem 
is in controlling the remainder in an asymptotic expansion of / in power of e. This 
requires that the macroscopic temperature gradient, i.e. |Ti — T2\/h, where h = eL 
is the thickness of the slab on the macroscopic scale, be small. 

Even if this apparently technical problem could be overcome we would still be 
left with the question of justifying the Boltzmann equation for such steady states 
and of course it would not tell us anything about dense fluids or crystals. In fact 
the Boltzmann equation itself is really closer to a macroscopic then to a microscopic 
description. It is obtained in a well defined kinetic scaling liniit. in which in addition 
to rescaling space and time the particle density goes to zeroO, i.e. A >> cr. 



3 Heat conduction in insulating crystals 

Excellent accountsL of the historical development of the theory of heat conduction 
in solids exist PtP^I so we will content ourselves here with some brief remarks. In 
(electrically) insulating solids, heat is transmitted through the vibrations of the 
lattice (in conductors the electronic contribution is in general much larger then 
the contribution due to the lattice vibrations). In order to use concepts of kinetic 
theory, it is useful to picture a solid as a gas of phonons which can store and 
transmit heat. In a perfectly harmonic crystal, the phonons behave like a gas of 
noninteracting particles and therefore the thermal current will not decrease with 
the length of the crystal placed between two thermal reservoirs. Thus a perfectly 
harmonic crystal has an infinite thermal conductivity: in the language of kinetic 
theory cr = and the mean free path A is infinite. A real crystal is not harmonic 
and, in the phonon picture, any thermal current will be degraded by the anharmonic 
forces in the lattice. Another source of finite thermal conductivity may be the 
lattice imperfections and impurities which will scatter the phonons and degrade 
the thermaLcurrent too. 

Debye lla devised a kind of kinetic theory for phonons in order to describe 
thermal conductivity. One assumes that a small gradient of temperature is imposed 
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and that the coUisions between phonons maintain local equilibrium. An elementary 
argument □ gives a thermal conductivity analogous to Eq. obtained in Section 
for gases, (remembering however that the density of phonon is itself a function of 
T) 

K ~ CyC^T . (7) 

Comparing Eq.(^ and Eq.(H) we see that p has been replaced by c„ the specific 
heat of the phonons, VT by c the (mean) velocity of the phonons, and A by ct, 
where r is the effective mean free time between phonon collisions. The thermal 
conductivity depends on the temperature yia r and a more refined theory is needed 
to account for this dependence. Peierls O used a Boltzmann type equation for 
phonons to investigate this problem. The Peierls theory singles out one phenomenon 
which gives rise to a finite thermal conductivityO. The momentum of phonons in 
collisions is conserved only modulo a vector of the reciprocal lattice. One can 
therefore classify the collisions of phonons into two classes: the ones where phonon 
momentum is conserved (the normal processes) and the ones where the initial and 
final momenta differ by a non-zero reciprocal lattice vector (the umklap processes) . 
Peierls theory may be summarized (very roughly) as follows: in the absence of 
umklap processes the mean free path and thus the thermal conductivity of an 
insulating solid is infinite. 

A success of Peierls theory is to describe correctly the temperature dependence 
of the thermal conductivity □. Furthermore, on the basis of this theory, one does 
not expect a finite thermal conductivity in f-dimensional mono-atomic lattices with 
pair interactions: this seems so far a correct prediction, see Section |l0|. _ 

The justification of the Boltzmann equation for phonons has been questioned Cj. 
Various alternative mechanisms have been proposed which would give rise to a finite 
thermal conductivity, but it seems fair to say that, so far, no better theory jof heat 
conduction in insulators has been proposed. As Peierls himself puts it Il3: "It 
seems there is no problem in modern physics for which there are on record as many 
false starts, and as many theories which overlook some essential feature, as in the 
problem of the thermal conductivity of [electrically] non-conducting crystals" . 

To find a mathematical description of thermal conduction in crystals we need 
to specify the Hamiltonian of the system or at least some appropriately idealized 
version of it. A model crystal is characterized by the fact that all atoms oscillate 
around given equilibrium positions. The equilibrium positions can be thought of as 
the points of a regular lattice in M . For simplicity we will assume that the lattice 
is simply Z . Although c? = 3 is the physical situation one can be interested also in 
the case d — 1,2. (The d — I system may show finite thermal conductivity without 
violating the Peierls criteria if we admit one particle, non momentum conserving 
interactions.) 

Let A C Z"* be a finite set and denote by N its cardinality. Each atom is 
identified by its position Xi = i + qi where i e A is the equilibrium position and 
qi (z R'^ is the displacement of the particle at lattice site i from this equilibrium 
position, and we denote by pi its momentum and m its mass. Since inter atomic 
forces in real solids have short range, it is reasonable to assume that the atoms 
interact only with their nearest neighbors via a potential that depends only on the 
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relative distance with respect to the equilibrium distance. 

As already noted it is useful to allow an external confining f-body potential 
which breaks the translation invariance. Accordingly the Hamiltonians that we 
consider have the general form 

^(^'^?) = E||+ E nqi-qj)+Ef^i(qi)-E|^ + ^(^?)' (8) 

ieA. |i-j| = l i ieA 

where P = (pi)igA and analogously for Q. We shall further assume that as |q| oo 
so do Ui{q) and V{q). The addition of Ui{q) pins down the crystal and ensures 
that exp [—f3H{P, Q)] is integrable with respect to dPdQ and thus the corresponding 
Gibbs measure is well defined. Observe that for many purpose it is enough tn nut 
the potential Ui on only some of the atoms, e.g. the ones on the boundary of ALjO. 
We note finally that when A C one can still consider that pi and qi € M"^, 
V ^ d, but wc will generally assume that S R''. 

Remark: While the Hamiltonian in Eq.(H) looks similar to that in Eq.(^) the 
meaning and domain of the Q variables is entirely different. In a fluid all the 
particle are identical and the particle with label i interacts with any other particle 
whose position q^ is close to q^, q/ G A C M'', I — \^ . . . . The pair interaction 
potential (/'(q) is of finite range, e.g. hard balls, or decays rapidly with distance. For 
the crystal in Eq.(||) q^ is the deviation from an equilibrium position i e A C Z'', 
etc. 



4 Microscopic models of heat reservoirs 



To produce a stationary heat flow in a system, be it a gas or a crystal, the system 
must be coupled to at least two heat reservoirs at different temperatures. A physical 
coupling is one which acts only at the boundary of the system leaving the dynamics 
in the bulk purely Hamiltonian. 

Since a realistic description of heat reservoirs and coupling is out of the ques- 
tion various model reservoirs have been used in analytical and numerical studies. 
We give here some examples which will be used later (other choices are of course 
possible). The expectation is of course that the different models will give the same 
behavior away from the boundary when the system is macroscopic. This has not 
been proven in any example, see p2| , ^ 



. 1 Stochastic 



reservoirs. 



We have already discussed one such model of reservoirs commonly used for fluids 
in Section ^. This corresponds to Maxwell boundary condition discussed in Eq.(||) 
for a gas in a rectangular slab. More generally given a fluid in a domain A C M'' 
a particle hitting the wall of the container confining the system at a point r e dA 
will bounce back into A with a Maxwellian distribution of momenta 

/r(dp) = ^P • n(r)e-'3(■-)^^dp , (9) 



luclO: submitted to World Scientific on February 5, 2008 



where h{r) is the inward directed unit vector normal to dA at r and /3~^(r) is the 
preassigned temperature at r. 

For sohds, which are usually not confined to any fixed spatial region by external 
walls, it is sometimes mathematically convenient to use Langevin type reservoirs 
which act on the atoms at the "edge" of the crystal. For definiteness we will choose 
A to be a chain of particles or a parallelepiped in higher dimension (with suitable 
boundary conditions), A = {i G Z ; I < ik < Nk , I < k < d}. We assume that the 
particles at the "left" boundary {ieA;ii = l} are coupled to a heat reservoir at 
temperature and that the particles at the "right " boundary {i e A; ii = A^i} 
are coupled to a heat reservoir at temperature Tr. We set C = Ni the length of 
the crystal and A = N2 - ■ ■ N^y its cross section. 

For the particles at the boundary of the crystal in contact with a heat reservoir, 
the Hamiltonian equations of motion arc modified by the addition of an Ornstein- 
Uehlenbeck process: 

TOiPi = -Vq,V(q) - A„pi/mi + (2A„T,)i/2^,(t) . (10) 

In Eq.(p^), a G {L, R} are indices of the reservoirs, Aq describes the strength of the 
coupling to the reservoir with temperature of the reservoir, and i,a{t) is a white 
noise, i.e. , a Gaussian random processes with covariance {S,a{i)(,i3{s)) = Saf3S{t — s). 
The form of the coefficients is chosen so that the dynamics satisfy detailed balance. 
This implies in particular that if the system is coupled to a single reservoir at 
temperature T, then the Gibbs measure with density cxp {-T'^H^P, Q)) IS a 
stationary state of the system. 

With any such a choice of stochastic process to model the reservoirs the dynam- 
ics is described by a stationary Markov process in the phase space of the system 



4-2 Hamiltonian reservoirs. 

In this case the reservoirs themselves are modeled by infinite Hamiltonian systems 
and the full system consisting of reservoirs+system is Hamiltonian^. An alternative 
but equivalent point of view is to start with an infinite system and to consider a 
finite subsystem of it as the system and the remaining part as the reservoirs. A 
non-equilibrium situation is obtained by choosing suitable initial conditions for the 
part of the total system which describes the reservoirs, e.g. the initial conditions 
of the reservoirs are assumed to be distributed according to a Gibbs measure with 
corresponding temperatures]^. 

Studied in ||77| , |80| , |69|] the simplest version of such a total system consists of an 
infinite chain. The "left" reservoir consists of the particles with labels in (—00, —N] 
and the right reservoir consists of the particles with labels in [N, +00). The system 



" A series of results has been obtained for systems coupled to a single reservoir, sunjpnias non- 
relativistic atoms coupled to the quantized electromagnetic |gi;14||ati temperature zero ErCrcl, finite 
level atoms coupled to a boson fiel|d-|at|positive temperature rT rir H riclassical particles coupled to 
a scalar field at temperature zero r'lP^ or at positive temperature Ej- 

''For quantum spin systems axioms are formulated in |ti'8| which establish the existence of a 
stationary state and its mixing property for finite systems coupled to several reservoirs at different 
temperatures 
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consists of the particles in the middle. At time t — Q, the reservoirs are assumed 
to be in thermal equilibrium at temperatures and Tr. This is clearly readily 
generalized to higher dimensions. 



4-3 Hamiltonian/ Stochastic reservoir. 

While Hamiltonian reservoirs are in principle the right ones to use they are totally 
intractable without further simplifications. When this is done it is actually possible 
to find examples in which one starts with a Hamiltonian reservoir and by integrating 
out over the degrees of freedom of the reservoirs, ends un with a stochastic evolution. 
Many models of this type have been constructecOEjO. We describe here a model 
considered in p5| , pT]j2^j20|] . The system is a finite chain of anharmonic oscillators 
coupled at each end to a reservoir modeled by a linear d-dimensional wave equation, 
which is the continuum limit of a d-dimensional lattice of harmonic oscillators. The 
dynamics of the infinite system, crystal+reservoirs, is Hamiltonian. One makes the 
statistical assumption that, at time t = 0, the reservoirs are in thermal equilibrium 
at temperatures Tl and Tr. Since the reservoirs are linear, one may integrate 
them out, and, by our assumptions on the initial conditions of the reservoirs, the 
resulting dynamics for the crystal is stochastic, though in general not Markovian. 
Nevertheless, the fact that the reservoirs are described by a wave equation, together 
with special choices of the coupling between the reservoirs and the chain, permits cil 
enlarging the phase space of the crystal with a finite number of auxiliary variables 
so that the dynamics is Markovian on the enlarged phase space. 

In the simplest case of coupling one variable per reservoir is enough and the 
resulting equations for the A'' oscillators are 

qi = -Vq,V(Q) + ri, 

- -Vq,V(Q), j = 2,...,7V-l, 
qn = -Vq„V(Q) +rfl, 

h = -M^R - >^W) + {2-iR\lTR)^/^WR{t) , (11) 

In Eqs.(|rT|) Al and \r describe the coupling strength to the reservoirs, 7^, 7_r are 
parameters describing the coupling and wl, wr are white noises. 

If the temperatures of both reservoirs are the same, = Tr = T, then the 
stationary state is given by the generalized Gibbs measure with density 

Z-'exp(^-^G{P,Q,R)^ (12) 

where Z is a normalization constant and the generalized "Hamiltonian" G is given 
by 



G(P, Q, R) = - qirL j + - qNTRj + H{P, Q) . (13) 

If one integrates the generalized Gibbs state, Eq.(p^), over the auxiliary variables 
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tl and rji one finds 

J dvLdvR Z-1 exp (-^G{P, Q, R)^ = Z'^ cxp (^-^H,ff{P, Q)^ , (14) 

where Hcs{P, Q) = H{P, Q) + A|qi/2 + A^q„/2 and Z a normalization constant. 
In view of this it is natural to consider Hcs{P, Q) as the energy of the chain. 

4-4 Thermostats. 

A fourth-. way of modeling the reservoirs is by deterministic (non-Hamiltonian) 
forcesLllEil. Such models of reservoirs are usually called thermostats. An exam- 
ple of such reservoirs which are widely used in numerical work, are the so called 
Nose-HooveiHEJ thermostats. Imposing these thermostats on small parts of the 
system^ (on the left and on the right) A^ and An, the equations of motion of 
particles in those region of the box are respectively 

mqi = -VqjV((5) - ClQi 

mqii = -Vq,V(g) - CflQi (15) 

where g Aq for a fluid and i G Aq. for a crystal, a G {L,R}. The variable (a 
model the action of the thermostat and satisfy the equations 

In Eq.(p^), 8 is interpreted as the response time of the reservoir and Ta is the 
temperature of the a-th reservoir. 

A limiting case of Eqs.(|l^)(|l6|) is when we let 8^0. This limit can be formally 
taken and the model becomes equivalent to the so called Gaussian thermostat. This 
means that one computes Cq, as a function of P and Q in such a way that the kinetic 
energy of the particles in A^ or A/} is a constant of the motion. After a simple 
calculation one gets for the chain: 

aiQ,p) ^ ^-<-- ' ^'^^ ' '^T' ' ^'^^ + , (17) 

I^KIl Pi 

and similarly for Cr{QiP)- Here f/(q) = — V?7i(q). One may also prescribe Gaus- 
sian thermostat in which the total energies, instead of just the kinetic energies, in 
Ai and A/j are kept fixed. 

5 Existence and Nature of Heat Conducting SNS 

Suppose we are given a system described by a Hamiltonian of the form Eq.(^) or 
Eq.(|8|) and that we have chosen a given model of heat reservoirs. We shall now 

'^It is also possible to give "mechanized" models of heat transport in which one imposes the 
presence of a heat flow through the application of a particular force on the bulk of the system. 
This strong modification of the Hamiltonian character of the dynamics seems unnatural to us, 
»ugh it can be useful for numerical simulations, and we will not discuss these models here 
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formulate a sequence of statements (of increasing mathematical difficulty) on the 
properties of the resulting dynamical system. 



5.1 Existence, Uniqueness and approach to the Stationary State. 

The first property that we want to prove is existence and if possible also uniqueness 
of a stationary state. For the case when all reservoirs are at the same temperature 
T existence is generally obvious - after all the reservoirs are chosen so that they 
leave the canonical Gibbs distribution or some variation of it invariant under the 
time evolution. Uniqueness and approach to this equilibrium state presents more 
of a problem and may not even be true for certain type of reservoirs and initial 
states. 

The real problem of interest for us is when the reservoirs are at different tem- 
peratures. We expect that if the dynamics is stochastic (e.g. for models 4.1 and 4.3 
of reservoirs) then "almost any" initial distribution of the state of the system con- 
verges to a unique stationary state which is mixing. This is however, in general, a 
mathematically non-trivial problem. The isolated system has a non-compact phase 
space and has many invariant states. The coupling to the reservoirs induces a drift 
towards a state determined by the reservoirs. Since however the coupling to the 
reservoirs occurs only at the boundary, the proof of the existence of an invariant 
measure requires a good understanding of how energy is transmitted through the 
system. There are in fact only few examples (to be discussed in Section |^) where 
this behavior has been proven. 

For general Hamiltonian or thermostated reservoirs the problem seems to be 
mathematically out of reach at the present time. Starting in a state that corre- 
sponds to the product of two equilibrium states for the reservoirs times a generic 
initial distribution for the crystal we then expect, that in the long time limit the 
marginal distribution for the system will approach some limit. When the two in- 
finite reservoirs are initially at different temperatures, the limiting state should 
describe a system having a temperature gradient and a heat flow. Observe that in 
general the state for each reservoir at times t > will not be the invariant state at a 
given temperature which is stationary for the isolated reservoir. It is this fact which 
makes the problem of general Hamiltonian reservoirs much more difficult that than 
of stochastic reservoirs. It is only in very special cases (essentially no interaction 



inside the reservoir) such as that discussed in subsection 4.3 where this can be dealt 
with. 

For thermostated systems the temperature of the reservoirs is already given by 
the equations of motion so we expect again to have a unique invariant distribution. 
Moreover these systems have the property that the phase space volume is not con- 
served by the dynamics so that, in general, no invariant measure will be absolutely 
continuous with respect to Lebesgue measure. In this case we need a criterion to 
chose the "physical" invariant distribution. A natural choice are the so called SRB 
states. These can be characterized by assuming that the system was in equilibrium 
in the very distant past and that at some point a forcing was switched on and this 
drove the system to a steady state distribution. More mathematically this means 
that we consider the weak limit of a probability distribution absolutely continuous 
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with respect to Lebesgue, e.g. the canonical distribution that characterize the sys- 
tem when all thermostats have the same temperature, under the time evolution, i.e. 
if is the flow describing the evolution of the thcrmostated system then a state /i 
is called SRB if it has the property that 

t^SRBidX) = lim $*A(dX) , X = {P,Q) (18) 

t — >oo 

where $* indicates the adjoint and X{dX) is the given initial distribution. We 
observe that we often cannot choose directly the Lebesgue distribution because 
the phase space of our system in not compact like for the Nose-Hoover thermostat 



Eq.(16) 



Although the definition Eq.([l8[) is interesting for its similarity to the ones used 
in the previous system-|-reservoirs models another characterization of these mea- 
sure is obtained by saying that they represent the statistics of the motion. More 
precisely, given any observable, the average of this observable with respect to the 
SRB distribution is equal to its time average along a trajectory starting from almost 
every point (with respect to Lebesgue measure). In formulae we can say that if ^ 
is the SRB distribution then 

1 /•* 

hm - / dtS^tix) = nsRB{dX) (19) 

t^OC t Jq 

for a set of X of full (or at least positive) Lebesgue measure. In Eq.(|l^) the limit 
is to be understood as a weak limitj^. 

5.2 Heat Flow and Entropy Production in Reservoirs. 

Since our interest here is specifically in Fourier's law our next question about the 
stationary state of the system coupled to two reservoir at different temperatures is 
the existence and nature of the heat flux across the system. It is clear that even 
existence is not automatic: most trivially just imagine that our system is composed 
of two noninteracting parts each coupled to a single reservoir. 

To study the heat flux through the system we first define a local energy density. 
For a fluid let q e A C M'' then 



N 



Mq;P,Q)=^<5(q-q,) 

i=l 



(20) 



2m ^-^ 

where the square bracket is identical to that in Eq.(|4|). For a crystal with nearest 
neighbor interactions we deflne the local energy density at site i as 

h{i; P, Q) = |L + ;7(qO + (^(^^ " ^-i^) + ^(^'+1^ " ^0) (21) 

k 

where Ifc is the d-dimensional vector with all components except the k-th equal 
to 1. 



''The only thermostated "physical" model for which a SNS corresponding to an SRB measure with 
the desired properties has been proven is the Moran-Hoover model of a single particle moving 
among fixed periodia-spatterers (see Section nh to which is added an external electric field and a 
Gaussian thermostato. 
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Given the local energy density we can define a local microscopic heat flow "if 
through the continuity equation. To avoid repetition we shall do so only for the 
crystal. Writing 

dh{i; P, Q) 



dt 



= V*(i) (22) 



where V*(i) = Efc<9».^'fc(i) with a,,*j(i) - (^^(i + Ifc) - «'j(i))/2. It easy to 
verify that 

M^.(i)^f(qi-qi+ijPi±|i±i^ (23) 

where f(q) = — Vy(q). We will usually be interested in the heat flow $(j') through 
the plane {i e A ; ii = j}. It is clear that we can integrate eg. (p^) and obtain 

■=(3.»2 id)' 

l<n<«i 

Of course the heat current inside the system in the steady state is just the 
energy flux from one reservoir to the other, presumably from the one with the 
higher temperature to the one with the lower one. Let S"* be the time evolution for 
observables (averaged over the realizations if the evolution is stochastic) . Since the 
equations of motion are Hamiltonian except at the boundary, one finds that the 
time derivative of the energy is given by 

^S'H{P, Q) = -5* ($L + $fl) (25) 
dt 

where <i>L {^r) depends only on the variables of the left (right) boundary of the 
system. It is natural to interpret $l as the flow of energy from the system to the 
left heat reservoir and similarly for We suppose that a stationary state /i exists, 
and, for any observable /, we set /i(/) = / fd^,. One obtains 

-^l{^^L+ $7?) = M {j^S*H{p, q)^ - , (26) 

and therefore 

^i{^L) = -^^{'^R) ■ (27) 

To check that the heat flux is iiideed as expected it is useful to define the entropy 
production a of the reservoirs asQ 

a=i^ + ^, (28) 

i.e., (7 is the sum of the energy flows into the reservoirs divided by the temperatures 
of the reservoirs. This (microscopic) definition of the reservoirs entropy production 
is in accordance with our notion of heat reservoir at specified temperature. It does 
not require that the system itself be close to equilibrium. A convenient way of 
proving that if Tl > Tji then heat is fiowing through the system from left to right 
is to show that in the steady state 
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and 



fi{a) =0 if and only if Tl=Tr. (30) 

(One may also consider the heat flow <l>(i) (or the corresponding J(q) for a fluid) 
inside the system and formally define a corresponding entropy production cr^ = 
(T^^ — Tj^^)<^{i). One obviously has that /z($l) = ^(<i>i) in the stationaxy^taie 
but (Ti is not the macroscopic entropy production density inside the systernEZnljIlj.) 

If Tl and Tr are close, one expects linear response theory to be valid. Setting 
T ={Tl+ Tii)/2 and ST = {Tl - T^), formal perturbation theory gives 

/•OC 

M(f) = / dtfiQ{aSQa)+ lower orders in , (31) 
Jo 

where fio is the equilibrium Gibbs state at temperature T and Sq the time evolution 
for observables with — Tji = T. For the heat flux $ — $(j') one obtains 

^($) = ^ ^ / dt^io{<^Sl^)+ lower orders in (5r . (32) 

It is important to note that in Eqs.(^ll) and (|3^), the reservoirs are still present via 
the time evolution 5*0. Although very similar this is not the Green-Kubo formula 
which will be discussed below. 

5.3 Fourier's Law. 

Assuming that 1. and 2. have been proved, we can then define the heat conductivity 
Kc as in Section ^ where C the length of the system (fluid or crystal) in microscopic 
units and A the area of its crossection. Since 5T/C is the average temperature 
gradient the heat conductivity at temperature T should then be given by 

- hm £ hm -1; (/i($)M) , (33) 

£— >oo aT— >0 01 

i.e., n — k{T) is the heat flux per unit area divided by the temperature gradient. 
One might also have taken the limit A oo in Eq.(^3|). As might be expected it 
is the limit £ — *■ oo which is the crux of the matter. 

6 Summary of Exact Results 

We now summarize briefly the limited number of results relating to points 1-3 of 
the last section. 

6.1 Fluid systems. 

Consider a system of N particle in A C M"^, d >2, with Hamiltonian given in Eq.(||), 
such that u(qi) = and the pair potential 0(|q|) is positive, with 0(0) = Ci, and 
-C2 < 0'(|q|) < 0, < Ci,C2 < 00, for \q\ > 0. Then it was shown in |3||],|| 
that, using Maxwell boundary conditions Eq.(^ with temperature T(r) > 0, r e dA 
(A a regular domain) there exists a unique stationary fi. Furthermore this /i is 
absolutely continuous with respect to Lebesgue measure on n = X R and is 
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approached, as t ^ oo, from almost any initial (P, Q), i.e. the set 57' for which the 
approach may fail has Lebesgue measure zero. 

The argument makes use of the boundedness of the force acting on any particle 
in the interior of A. This assures that any particle with a sufficiently high speed 
will hit the boundary with only little deviation from a straight path. This and the 
fact that the force — 0'(|q|) is everywhere positive insures effective contact between 
the system and the stochastic boundaries which, according to Eq. ^ , "spread" the 
velocity of particles which hit it. This yields something like a "Harris condition" 
which guarantees existence, uniqueness and approach to the stationary state. 

Using general technique developed in Eq.(|2^) for systems in contact with 
stochastic reservoirs satisfying detailed balance is immediate. It is probably also 
possible to prove inequality ( |30| ) for such systems but the latter has not been done 
generally as far as for we know; see below. 



6.2 Harmonic Crystal. 

A system with Hamiltonian given by Eq.(H) in which both V and Ui (when it does 
not vanish) are quadratic functions of their arguments is an ideal harmonic crystal. 
When such a system is placed in contact with stochastic reservoirs of the Langevin 
type the resulting process and thus also the stationary measure is Gaussian and 
one only needs to compute the covariances. This was done essentially explicitly for 
a chain in [75| ]. The most important difference with the equilibrium state is that 
there are now non- vanishing covariances between position and momentum variables 
proportional to ST. One finds uniqueness and approach to the stationary measure 
fi that satisfies Eqs.jp^) and ( |3C| ) 

As already mentioned however the heat flux /i(<l>) is essentially independent of 
C and Kc defiiLcd_in Eq.(^3|) grows as £. (For the case of "random masses" kc 
grows as V31jO'L3). The solution of |Q was extended to d > 1 in j67| where there 
were also considered various possibilities for Ui, e.g. pinned down everywhere, only 
at the boundary or nowhere. Looking at the invariant measure in the limit C oo 
one finds that the decay of the position- momentum covariance is rapid when Ui ^ 
(at least on the boundary) but does not decay at all when Ui = 0. 

The case of an infinite harmonic chain with left and right portions acting as 
Hamiltonian reservoirs was investigated in The results are qualitatively the 

same as for the stochastic reservoirs: the heat current remain proportional to the 
initial temperature difference as t oo and the system approaches its stationary 
state which is again a Gaussian measureFI. 



6.3 Anharmonic Crystals. 

The anharmonic crystal coupled to Hamiltonian/Stochastic reservoirs described by 
Eq.(ll) has been investigated in |pT| , p2| , po| . JTechnical conditions on the growth a± 
infinity of the potential are needed: either til that V is quadratic at infinity or E3 
more general polynomial growth. (In the latter case, the one-body potential U 



■^For an infinite quantum harmonic chain with a special particle subject to a sufficiently small 
non harmonic potential, the exis tenca of stationarjj-states and their mixing property has been 
established, both for KMS states E3Ea and for SNS El. 
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grows more slowly at infinity then the two-body potential V). One assumes alaci 
that the two-body potential is strictly convex, and this condition alone implies E3 
that the stationary state is unique. Under these conditions the following results 
hold: 

1. Existence and uniqueness of the stationary state /j,. The stationary state is 
mixing, i.e. any initial distribution will converges to the stationary state as 
t —>■ oo. The stationary state has a C°° density which decays at infinity at 
least as fast as a Gibbs state with temperature equal to the maximum of the 
temperature of the reservoirs. 

2. The stationary state is conducting: One has = if and only if = 
and h{^r) > if Tl > Tr. Linear response theory is vahd: For a large class 
of observables /, the expectation value /^(/) is a real-analytic function of the 
temperature difference 6T. In particular, near equilibrium, one obtains, with 
T = (Ti + Tr)/2 

K^B.) = ^^D + 0{ST'), (34) 

Eq.^for the coefficient D has not been proved, but rather the slightly weaker 
form 

D = ^ioi<^>R.iLo^<^>R)), (35) 

where /iq is the Gibbs state Eq.(^) with temperature T and Lq is the genera- 
tor of the Markovian semi-group Sq associated with the stochastic differential 
equations (|ll|) with = Tr = T. Notice that, formally, one has 

dtSl (36) 

and, inserting Eq.(|3^) into Eq.(|35|) yields Eqs.(^). In order to prove Eq. (|36|) , 
one needs presumably some information on the decay of correlations. This has 
not been obtained so far. Nothing is known about the dependence of D on C 
and thus ipso facto about the validity of Fourier's law. 

7 The Green-Kubo Formula 

It is clear that, aside from the case of the harmonic crystal, which does not satisfy 
Fourier's law, none of the exact results quoted in the last section says anything 
about the local structure, e.g. about local equilibrium in the SNS. This means in 
particular that at this time we have no rigorous way of relating the local heat 
flux ^{^) to the gradient of the local temperature as defined in Eq.(^). Even in 
the absence of LTE one can define a local kinetic temperature by means of the 
average local kinetic energy. Thus for the crystal at site i, T(i) = ii{p? /Tn)/d. 
For the harmonic crystal T(i) in found to be uniform away from the ends, i.e. 
there is no temperature gradient. Even accepting such a definition of temperature 
(in numerical simulation, to be discussed later, T(i) is one of the most directly 
measured quantities) we are completely lacking at this point any rigorous or even 
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formal connection between the k defined in Eq.(p3|) and the usual Green-Kubo 
formula for the conductivity which is defined in terms of the time evolution of an 
isolated system in equilibrium. Denoting by Sq the HamiLtniuan evolution of the 
isolated system the thermal conductivity kgk is given by e3E3 

where (•) denotes the microcanonical average and the energy density is chosen such 
that it corresponds to the thermodynamic energy at temperature T (since we are 
in equilibrium this is the same as the kinetic temperature). If the total momentum 
n is conserved it has to be set equal to zero. Alternatively one may use in Eq.(^7|) 
truncated correlation functions (*i>S'o''')l;,n — {^^o^)E.n — {^)E.n and then average 
over E and 11 using the canonical distribution^. One expects that the equivalence 
of equilibrium ensembles will extend also to this case. 

The Green-Kubo formula, Eq.(^), also makes sense for a system with a few 
degree of freedom where it can be related to the vpjiance in the fluctuation of the 
current; see the article of Bunimovich and Spohn O for a discussion. There is no 
clear connection however between the integral in Eq.(p7|) for a small system and 
the K in Fourier's law. From a mathematical point of view it is not even clear how 
to prove equivalence for macroscopic systems, i.e. show that k = kgk- It would be 
nice to find even a formal argument establishing the equivalence. 



8 Entropy Production and Large Deviations 

The proper definition of microscopic entropy production has attracted much at- 
tention in recent years. The interest comes from the observation of an interesting 
symmetry property in the large deviation functional associated to the phase space 
volume contraction rate in thermostated systems. This property was first observed 
numerically in and then proved under strong hyperbolicity condition in |]3l|] . 
For such systems the phase space volume contraction has strong connection to the 
entropy production in Section |^. Using this connection the fluctuation theorem has 
been e xtend ed to large deviations of the entropy production of various stochastic 
systema£3'E3. For crystals with stochastic reservoirs or for the model considered 
in Section 4.3, only formal proofs of the Gallavotti-Cohen fluctuation theorem are 
available so far. 

The Gallavotti-Cohen fluctuation theorem can be formulated as follows: In 
both deterministic and stochastic systems one identifies an observable cr{X) as the 
"entropy production" although its identification with the thermodynamic entropy 
production is not always clear. If one now considers the ergodic average, 

1 /•* 

at (X)^- dsa {X{s)) (38) 

^ JQ 

where for deterministic dynamics X{t) is trajectory of the system while for stochas- 
tic dynamics X{t) is a particular realization of the random process. Then by the 

/If one does not fix 11 in an equiUbriup*] ensemble for which the total momentum is conserved 
then the integral in Eq. (BTI) is divergento but this does not say anything about hqk- 
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large deviation principle (assumed to hold) there exists an e{p) such that for any 
interval / 



lim logProb(crt(x) £ I) = inf e{p) . (39) 

t^oo pel 

The fluctuation theorem asserts that the odd part of e{p) is linear with slope —1, 
i.e. e(p) — e{—p) = —p. 

In the deterministic case the entropy production variable that one consider is, as 
already mentioned, the phase space contraction rate which, in the case of Gaussian 



thermostats, coincides with entropy production as we defined it in Subsection 5.2 



This can be easily seen by computing the divergence of the equations of motion 



quoted at the end of Subsection 4.4. Already for the case of the Nose-Hoover ther- 
mostat the phase space contraction and the entropy production are not the same 
quantities although their average value in the SNS are the same. For anharmonic 
chains with a two body potential given V{x) = + /3x* the validity of the fluctua- 



tion theorem has been checked numerically in 1 59 1 for the entropy production (but 
not for the phase space contraction). 

The connection between the Gallavotti-Cohen fluctuation theorem and its var- 
ious generalizations can be related to the following observation: For a stochastic 
model described by a Markov process one can consider the measure P on the path 
space induced by the evolution. Let be a path leading from X to F in phase space 
in the time interval [0, t] and consider now the transformation which maps the path 
(j) into the time reversed path leading from lY to IX, where / is the involution on 
phase space which reverses the velocities of all particles. This transformation maps 
the measure P into a new measure P which is absolutely continuous with respect 
to P with Radon-Nykodym derivative given by 



^(</)) = exp (^Rim) - Rim) + ^ <^iHs))ds^ , (40) 

where a is the entropy production defined in Eq. (p8|) . Eq. ( ^0|) states, roughly 
speaking, that the probability of a time reversed path is equal (up to a boundary 
term) to the probability of taht path times the exponential of the integrated entropy 
production along thiSpflath. This property can be assumed as a general definition 
of entropy productions and can be seen as a generalization of detailed balance. In 



fact in the equilibrium case (Tl = T/j = T), the right hand side of Eq.(40) depends 
on (f) only through its endpoints and is equal to exp (T^^(i?(0(t)) — H{(f){Q)))) and 
this is precisely detailed balance. 

For thermostated systems the above calculation is quite delicate because, al- 
though one can still consider the measure P on the path space, it will not be 
absolutely continuous with respect to the time reversed measure P. The original 
proof of the fluctuation theorem in can be thought of as based on the construc- 
tion of a series of approximants to the measure P which satisfies a relation formally 
similar to Eq.(En). 
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9 Local equilibrium 

Since we have nothing to say about LTE in Hamiltonian heat conducting systems 
we discuss briefly some models where the Hamiltonian dynamics in the bulk is 
modified by the addition of stochastic forces. These models show that a mechanism 
sufficiently strong to destroy the coherence of the phonons does indeed produce both 
a temperature gradient and a normal thermal conductivity, even in one dimensional 
systems. _|_. 

The first such model EHIa is a harmonic chain coupled to self-consistent reservoirs. 
Each particle of the chain is coupled to its own stochastic reservoir at temperature 
Ti, i = 1, ■ ■ ■ , N , modeled by Langevin dynamics. The reservoirs coupled to the par- 
ticles in the bulk supposedly simulate the effect of strong anharmonic interactions. 
The temperature of the reservoirs coupled to the first and last particles are at fixed 
temperature Ti — Tj^ and T^v = Tr. The temperature of the remaining reservoirs is 
fixed by a self-consistency condition: one requires, that in the stationary state there 
is no net epergy exchange between the reservoirs and the particles in the bulk. It 
is argued El that in the limit of largeJV, the system exhibits a temperature gradient 
and Fourier's law i^ndeed satisfiecO. 

Another model ED is a one-dimensional chain of quantum mechanical atoms (each 
of which has a finite number of energy levels). The first and the last atoms of the 
chain are each coupled to heat reservoirs at temperatures Tj, and Tr. The atoms in 
the chains are not directly coupled to each other, but each pair of nearest neighbors 
in the chain is coupled to a heat reservoir. Since each atom has only a finite number 
of energy levels, it is possible to choose the coupling to the intermediate reservoirs 
in such a way that no energy is exchanged between the particles in the bulk and 
the reservoirs. Compared to the previous model, no self-consistency condition is 
needed and the transfer of energy between the intermediate reservoirs and the 
system always vanishes, not only in the stationary state. The model is studied 
in the weak coupling (or Van Hove) limit, where the coupling between the chains 
and all reservoirs goes to zero and the time is suitably rescaled. In this limit the 
evolution is described by a quantum semigroup. For a suitable choice of couplings 
the model can be solved exactly and it exhibits a temperature gradient and Fourier's 
law. _ 

Another model considered c3 consists of a chain of uncoupled harmonic oscilla- 
tors {Ui{qi) = <Z,f/2, V{qi — Qj) = 0). The oscillators at the boundary are coupled 
to heat reservoirs modeled by Glauber processes which thermalize the oscillators 
according to the Gibbs distribution at temperatures Tl and Tr. The energy is 
exchanged between the oscillators in the chain according to the following (micro- 
canonical procedure): at each pair of nearest neighbor sites, there is a clock with 
exponential law, when it rings the energy of the pair of particle is redistributed 
in a uniform way, keeping the total energy of the system constant. This model is 
exactly solvable and satisfies Fourier's law. 
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10 Numerical results 

The availability of fast computers has permitted the investigation of heat current 
carrying SNS via numerical simulations. Many such simulations have been per- 
formed, see e.g. | |T^ , |73| , ^ , |60| , ^2|j4C| Jl|,|3^ for recent works on various models. Earlier 
works were not always consistent but recently jCareful simulations of both one and 
two dimensional crystals have beccpae availableE3. A coherent picture now seems to 
emerge from the numerical results^. We describe this briefly and refer the reader 
to | ]6]| for a more complete overview. 

For obvious reasons, it is convenient in numerical works to keep the number of 
degrees of freedom as low as possible. On the other hand stochastic differential 
equations, like the ones used in the Langevin stochastic reservoirs require very 
good random number generators and present numerical problems connected with 
the singularity of the covariance of the white noise. For this reason the simulations 
done using the deterministic thermostat are easier and likely more reliable. In 
any case there is good agreement (in the data if not always in the interpretation) 
between simulations with different thermostats. 

Since it is not easy to look for invariant measures with numerical simulations 
one typically computes the time average of a few interesting observables along a 
given trajectory. Assuming the validity of Eq.([l9|) this represents averages with 
respect to the invariant measure. This permits investigation of questions like the 
temperature profile in a chain in the steady state or the value of the conductivity 
for which, as we have seen, there is no analytical result at this time. We will focus 
mainly on this last question. 

1. In one dimension, where most of the simulations have been conducted, the 
conductivity, when f/i(q) = 0, appears to behave, as a function of the length 
£ of the chain, as with a a positive exponent, a = 0.4 for the anharmonic 
chain. On the other hand if C/i(q) ^ (typically one consider C/i(q) = ^uj'^q^) 
one finds a finite conductivity if some nonlinearity is present in the system 
(for the linear case a = lEj). In this situation the exact form of the non 
linearity seems irrelevant, analogous result are found adding a 4-th order term 
to C/(q) as well as to V{q). Similar results are obtained if one compactifies 
the configuration space of Q, e.g. by considering each q as a point on a torus, 
obtaining what can be called a chain of rotators. In all cases the system has 
a well defined, approximately linear, temperature profile, although there can 
be a finite jump between the temperature of the first and last oscillator and 
the temperature of the respective thermostats. Moreover for one dimensional 
systems simulations using different kinds of reservoirs appear to yield similar 
values of the exponent a. 

2. Recent simulations seem to show that the conductivity in two dimensions is 
logarithmically divergent if J7i(q) = 0. Although it is not easy to see a loga- 
rithm is such a situation the simulation in strongly suggests this conclusion. 
Moreover in this case one can try to compute the thermal conductivity hqk as 
given in Eq.(^7|)p| and compare it with the n obtained from Eq.(|3^) (without 

®In this case the equihbrium ensemble used is {■)e,-w described in the comments after Eq.(p7|) 
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the limit £ — > oo). The agreement obtained from the numerical data support 
the validity of the equality k = kqk- Although we do not know any results 
on this direction we believe that adding a confining U{q} to this system will 
make the conductivity finite. 

3. We further expect that the conductivity in three dimensions will be finite, with 
or without any on site potential. 

The above picture can be interpreted in term of the Peierls theory which, as 
mentioned in Section ^ relies on umklap processes to produce a finite conductivity. 
We also note that in one and two dimension, there is no stable crystal without an 
on site potential, but there is localization in d > 3. The variance in the deviation 
of an atom near the center from its equilibrium position, when the boundary atoms 
are tied down, grows like C in d = 1, like log£ in rf = 2 and is finite in c? = 313. 
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